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Abstract 

A general formula for the canonical partition function for a system obeying 
any statistics based on the permutation group is derived. The formula ex- 
presses the canonical partition function in terms of sums of Schur functions. 
The only hitherto known result due to Suranyi [ Phys. Rev. Lett. 65, 2329 
(1990)] for parasystems of order two is shown to arise as a special case of our 
general formula. Our results also yield all the relevant information about the 
structure of the Fock spaces for parasystems. 
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There exist two approaches to parastatistics. The first is the field theoretic approach 1-3 
based on a generalization of the creation-annihilation operator algebra for bosons and 
fermions. This is the way in which parastatistics was first introduced by Green. 1 In the 
second approach pioneered by Messiah and Greenberg 4 and further investigated by Hartle 
Stolt and Taylor, 5 parastatistics arises in the quantum mechanical description of an assembly 
of iV-identical particles. Here the permutation group Sn plays a central role. Of the two, for 
calculational purposes, the first seems to have found greater favour with the workers in this 
field. 6 ' 7 Thus, for instance, using this approach Suranyi 7 has given the canonical partition 
function for a parabose and parafermi gas of order p — 2. The calculation involves a clever 
use of the simplifications which occur in the para algebra when the order of the parastatis- 
tics is two. In this letter we show that the quantum mechanical approach to parastatistics 
when combined with the machinery of symmetric functions 8 yields a powerful method which 
enables one to answer, with great facility, all questions pertaining to statistics based on the 
permutation group. In particular, we give the canonical partition function for a parabose or 
parafermi system of an arbitrary order. 

To set the notation we begin with a brief summary of some familiar results. This is 
also necessitated by the fact that it is an in-depth appreciation and examination of what is 
all too familiar along with the intution gained in the process which leads us, in one stroke 
to the desired results. Consider a Hilbert space H built by an iV-fold tensor product of a 
Hilbert space H of dim M. We shall assume that M > N. Let 1, 2, 3, • • • , M denote the 
basis vectors of H. The M N basis vectors of Ti correspond to each term in the product 

(1 + 2 + ■ ■ ■ + M)(l + 2 + ■ ■ ■ + M) ■ ■ ■ (1 + 2 + ■ ■ ■ + M) 

N factors ' ^ ' 

One may consider two decompositions of this set of M N states. 

1. Decomposition based on occupation numbers: 

Here one groups together states which have the same number of l's, 2's • • • etc., regardless 
of their location in the product. Each such group is characterized by a set of occupation 
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numbers which give the number of times 1,2, • • ■ M occur in the states in that group. All 
relevant aspects of this decomposition are encapsulated in the following decomposition of 
the symmetric function {x\ + • • • + %m) n , 

Z l ^ f (x 1 ,---,x M ) = (xi + x 2 Vx M ) N = J2 x , N \ , m x (x ir --,x M ) • (2) 

a Ai! • • • Am- 

\X\=N 

Here A = (Ai, A2, • • • , Am), Ai > A2 > A3 • • • > Am is a partition of A" (indicated by 
|A| = N) and m\(xi, ■ ■ ■ ,%) denotes the monomial symmetric function 8 corresponding to 
the partition A 

m x (x 1 ,---,x M ) = Y, x i 1x 2---x X m ■ (3) 

The sum on the R.H.S. of (3) is over all distinct permutations of (Ai,---,Am)- Setting 
x\ = xi = xm = 1 in (2) we obtain 

M Af = E x N \ ^ a(1,---,1) , (4) 

which tells us that each partition A = (Ai,---,Am) corresponds to m\(l, • • • ,1) sets of 
occupation numbers obtained by distinct permutations of Aj's and each such set contains 
AH/Ai! • • • Am' states. The number m A (l, •••,!) is given by 

m A (l,l,---,l)= - , (5) 

m 1 !m 2 ! • • • 

where m^s denote the number of times Aj's occur in the given partition A. 
Note also that with the identification 

x 1 = e _/3ei , x 2 = e"^ 2 , • • • , x M = e~^ M , (6) 

where ei, e 2 • • • , €m are taken to denote the "single particle energies" corresponding to the 
states 1, • • • , M, the symmetric function Z 1 ^ (xi, • • • , %) represents the partition function 
of the system under consideration. This fact is made more explicitly by rewriting (2) as 

Z™ f (x ir --,x M ) = (x l + --- + x M ) N = E -^T^ rX^xr---x n M M . (7) 

n^.nj. ■ ■ -n M ' 

(The symbol Z^ was introduced with a view to emphasizing this point. The superscript 
"inf" indicates that we are dealing with the infinite statistics 9 ) 



1. Decomposition based on the permutation group: 

In this decomposition we regard the M N states as the carrier space for an M N dimen- 
sional representation of the permutation group Sn- This reducible representation can be 
decomposed into the irreducible representations of Sn which, as is well known, are in one 
to one correspondence with the partitions of N. In this case, the relation analogous to (2) 
which summarises all relevant features of this decomposition is 

Z*^ f (xi,---,x M ) = (xi-\ \-x M ) N = n(\)Sx(xi,---,x M ) , (8) 

A 

|A|=Af 

where S\(xi, • ■ ■ , %) are the Schur functions 8 

det(x Xj+M ~ j} 

S x (x 1 ,---,x M )= I l<i,3<M , (9) 

det [Xj ) 

and n(A) denotes the dimension of the irreducible representation A of Sn- Setting x\ = x 2 = 
■ ■ ■ = xm = 1 in (8) we get 

M N = n(A)5 A (l,...,l) , (10) 

A 

|A|=JV 

which tells us that ^(1, • • • , 1) is the number of times the irreducible representation A occurs 
in this decomposition. 

So far we had been dealing with Ti. Following refs. 4 and 5 we now construct out of it 
a generalized ray space Ti p h y by 

• (a) admitting only those operators on Ti which are permutation symmetric, 

• (b) identifying those states in Ti which have the same expectation values for all per- 
mutation symmetric operators. 

These assumptions, together with the Schur's Lemma, imply that all states belonging to 
an irreducible representation A of Sn count as one state of Ti p h y - This immediately tells us 
that the symmetric function of degree N appropriate to Ti p h y is simply obtained by setting 
n(X) = 1 in (8). 
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(xi, - ■ ■ ,x M ) = S x {x ir ■ ■ ,X M ) . 



(11) 



A 



|A|=JV 



This is the key result of this work. (Here we use the superscripts HST to denote Hartle Stolt 
and Taylor in honour of their contributions to parastatistics). This symmetric function 
captures all aspects of TL p h y much the same way as (x± + • • ■ + xm) N does for TC. 

So far no restrictions have been put on A — the sum on the R.H.S. of (11) is over all 
partitions of N. We shall refer to this statistics as HST statistics. Parabose case of order p 
arises when we retrict the sum in (11) to only those partitions of N whose length /(A) (the 
number of the non-zero Aj's) is less than equal to < p. In terms of Young tableaux, this 
amounts to retaining only those irreducible representations pf Sn in which the number of 
boxes in the first column is < p. The appropriate symmetric function for this case is 



Similarly parafermi case of order p arises when we restrict A in (11) to those partitions whose 
conjugate partition A' is of length < p. In terms of Young tableaux this implies retaining 
only those irreducible representations of Sn in which the number of boxes in the first row 
is < p. The symmetric function appropriate to this case is 



Likewise, for the (p, q) statistics the corresponding symmetric function Z ( ^' q '(xi, • • • , xm) is 
obtained by restricting the sum in (11) to those partitions for which /(A) < p and /(A') < q. 
It may be noted that if p, q > N, all these cases reduce to HST. 

The symmetric functions above contain all the relevant information about the appropriate 
Hphy For instance the dimension of TL p h ys is obtained by setting x± — x-i — • • • xm = 
1. The appropriate canonical partition function is obtained by setting x\ = e _/3ei , x 2 = 
e~ l3e2 , ■ ■ ■ ,xm = e _/3eM . Some formulae which prove to be extremely useful in carrying out 
the sums in (12) and (13) with restrictions on the lengths of the partitions are as follows. 8 



Z N - B (x ir ■ ■ ,x M ;p) = Sx(x ir ■ ■ ,x M ) ■ 



(12) 



X \\\=N 
l(X)<p 




(13) 



5 



S x (x 1 ,---,x M ) =det(h Xl -i+j) ; l<i,j<l(\) , 
S\(xi,---,x M )=det(ex'.-i+j) ; l<i,j <l(\') . 



(14) 
(15) 



Here the complete symmetric functions h r (xi, ■ • ■ , %) and the elementary symmetric 
functions e r (xi, ■ ■ ■ , %) are defined as follows 



h r (x u • • • , x M ) = m A(^i, • • • , im) , 

A 

|A|=r 

e r (a^x , • • • , Xj[4 ) ^ ' 2-12 ' ' ' -Eir 



(16) 
(17) 

U<12<---<V 

Using these formulae one can express the appropriate symmetric functions in terms of either 
fa's or e's. As an illustration, let us consider the Bose case. Here, since /(A) < 1, we have 
only one term on the R.H.S. of (13) corresponding to A = (iV, 0, • • • , 0). Using (14) we 
obtain 



Z^(xi, ■ ■ ■ , xm) — h}sr(xi, • ■ • , xm) 



(18) 



Similarly, for the Fermi case, one has 



Zn(xi, ■ ■ ■ , xm) — &n(xi, • • • , xm) 



(19) 



Consider parabose of order 2. Using (12) and (14) we obtain 

Zx B (xi, ■ ■ ■ , x M ; 2) = h N + det 
which on simplification leads to 



A 1 +A 2 =iV 
A X >A 2 



y h\ 2 -i h\ 2 j 



(20) 



(xi, ■ ■ ■ ,x M ; 2) 



4(xi,---,x M ) if AT = 2P , 



h P+1 (x ir ■ ■ ,x M )h P (x 1 ,- ■ ■ ,x M ) if iV = 2P+l 



(21) 



The result for parafermi of order 2 is obtained by replacing fa's by e's. Thus we obtain the 
the results due to Suranyi which arise as a special case of (12) and (13). One can carry out 
similar calculations for any order p. 
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For the HST case a number of interesting formulae can be derived from the following 
result for the Schur functions 8 

E E s ^---^) = Uj^U7^ r) • ( 22 ) 

\\\=N 

The R.H.S. of (22), with appropriate identifications of the be seen to be the grand 

canonical partition function for the HST statistics. Further, setting x\ — • • ■ — xm = t and 
reading off the coefficient of t N in the resulting expression on the R.H.S. one obtains 



dim(H phys ) = ^(1, ••-,!) 



A 

|A|=JV 

[N/2] 

E 

s=o 



( M + N - 2S-1^ 
N-2S 



^ M(M - l)/2 + S - 1 ^ 

v s / 



(23) 



This formula is interesting in its own right as it gives the number of irreducible representa- 
tions of Sn which occur in the decomposition of the M N dimensional reducible representation 
of Sn discussed above. 

Finally, having derived the symmetric function appropriate to H p h y (with or without 
restrictions on the irreducible representations of the permutation group) we can immediately 
obtain all the relevant information regarding its occupation number decomposition. All that 
needs to be done is to expand the appropriate symmetric functions in terms of the monomial 
symmetric functions m,\(xi, • • • , %). Now since S\(xi, • • • , %) and m\{x\ ■ • • %); |A| = N, 
serve as bases for symmetric polynomials in xi, ■ ■ • , % of degree N, we can expand one in 
terms of the other 



S x (x ir ■ ■ ,x M ) = K x xm\(x!,- ■ ■ ,x M ) , (24) 

A 

|A|=JV 

where K x \ are Kostka-Foulkes numbers. 10 Substituting this in (11) we obtain 

Z§ ST (x u - ■ ■ ,x M ) = ^C^K xX )m x (x l} - ■ ■ } x M ) ■ (25) 

a x 

(For parabose or parafermi or (p, q) case one has to put appropriate restrictions on the 
partitions x m (26)). The coefficients Y, X K X \ immediately give us the number of states 

7 



corresponding to the occupation numbers A = (Xi,---,Xm) or any distinct permutation 
thereof. 

To conclude, by combining the approach propounded in refs. 4 and 5, with the theory of 
symmetric functions, we have been able to obtain partition functions for all statistics based 
on the permutation group. Detailed analyses of the thermodynamic properties derivable 
from these results would be published elsewhere. 
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